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Abstract 



Path integral methods are used to derive a general expression for the en- 
tropy of a black hole in a diffeomorphism invariant theory. The result, which 
depends on the variational derivative of the Lagrangian with respect to the 
Riemann tensor, agrees with the result obtained from Noether charge meth- 
ods by Iyer and Wald. The method used here is based on the direct expression 
of the density of states as a path integral (the microcanonical functional inte- 
gral). The analysis makes crucial use of the Hamiltonian form of the action. 
An algorithm for placing the action of a diffeomorphism invariant theory in 
Hamiltonian form is presented. Other path integral approaches to the deriva- 
tion of black hole entropy include the Hilbert action surface term method 
and the conical deficit angle method. The relationships between these path 
integral methods are presented. 
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I. INTRODUCTION 



Noether charge methods have led Iyer and Wald [|T],[2j to the discovery of two elegant 
expressions for the entropy of a stationary black hole in a diffeomorphism invariant theory 
in D spacetime dimensions. The first expression is 

S BH = 2vr / Q[t] , (1.1) 

where Ti. denotes the black hole bifurcation surface and Q[t] is the Noether charge (D — 2)- 
form associated with the horizon Killing field t a . The second expression is 

S BH = -2tt / d D - 2 x^e ab e cd Ut d , (1-2) 
Jn 

where a is the determinant of the metric on Ti, e a b is the binormal of Ti, and UQ bcd is 
the variational derivative jp of the Lagrangian with respect to the Riemann tensor TZ a bcd- 
Equation (1.2) is a generalization of the result obtained in Ref. || for black hole entropy in 
a theory described by a Lagrangian that depends on at most first derivatives of the Riemann 
tensor. The equivalence of expressions (1.1) and (1.2) is demonstrated in Ref. 0. 

More recently, Iyer and Wald |J and Nelson || compared the Noether charge approach 
with various path integral derivations of black hole entropy. The path integral methods 
all originate, ultimately, with the observation made by Gibbons and Hawking [0] that the 
partition function for the gravitational field can be expressed as a path integral.^] These path 
integral methods were developed within the context of specific theories, such as Einstein 
gravity or Lovelock gravity. They include (i) the direct expression of exp(S BH ) in terms of 
a path integral (the microcanonical functional integral) |p|,|T0|| ; (ii) the expression of S BH in 
terms of the Hilbert action surface term [0,0] > an d (hi) the derivation of S BH in terms of a 



nonclassical spacetime with a conical singularity |12|,|13[. Using the language and techniques 
of the Noether charge formalism, Iyer and Wald showed that the path integral methods 
(i) and (ii) yield the result (1.1), the black hole entropy expressed as the integral of the 
Noether charge Q[t], when applied to an arbitrary diffeomorphism invariant theory. Nelson 
has analyzed the relationship between the path integral method (iii) and the Noether charge 
result (1.1). 

Section 4 of this paper contains a direct derivation of the result (1.2), the black hole 
entropy expressed in terms of the variational derivative of the Lagrangian. This derivation is 
based on the microcanonical functional integral method (i) and bypasses the Noether charge 
formalism altogether. There may be some advantage to this. The Noether charge formalism 



The original calculation of Gibbons and Hawking was inconsistent. Their result for the partition 
function implies a negative value for the heat capacity. On the other hand, general arguments 
show that the heat capacity is necessarily positive for any system that can be characterized by a 
partition function. This problem was overcome by York H who showed that the partition function 
yields a positive value for the heat capacity if the boundary conditions in the path integral are 
imposed at a finite spatial location. 
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is a useful tool for deriving the first law of black hole mechanics but, by itself, it does not 
provide a logically complete derivation of black hole entropy. In order to extract the black 
hole entropy from the first law of black hole mechanics one must, in principle, supplement 



the Noether charge analysis with the quantum field theory scattering calculation [[lj] that 
leads to the identification of surface gravity (divided by 2tt) with black hole temperature. 
On the other hand, the path integral approach, although formal, does provide a logically 
complete framework in which black hole entropy can be derived and analyzed. For this 
reason, insights into the mysteries of black hole entropy, such as its statistical origin, should 
be obtained more easily from within the path integral formalism. 

The microcanonical functional integral method (i) is reviewed in Sec. 2. The derivation of 
Eq. (1.2) in Sec. 4 makes crucial use of the Hamiltonian form of the action, which is derived 
in Sec. 3. The three path integral methods mentioned above, having a common origin, 
are closely related to one another. The logical connections between the microcanonical 
functional integral method (i), the Hilbert action surface term method (ii), and the conical 
deficit angle method (iii) are discussed in Sec. 5. 

An important part of the present analysis is contained in Sec. 3, where an algorithm is 
developed that allows one to place the action for a diffeomorphism invariant theory in Hamil- 
tonian form. Specifically, it is shown that any action can be placed in "almost Hamiltonian" 
form, which differs from a true Hamiltonian form by the presence of extra undifferentiated 
variables (referred to as the x's) in the Hamiltonian constraint C±. A true Hamiltonian form 
of the action is obtained when the x's are eliminated through the solution of their algebraic 
equations of motion. If the rank of the matrix formed from the second derivatives of C± 
with respect to the x's is not maximal, then there are constraints on the canonical variables. 
These constraints must be added to the action via Lagrange multipliers. In practice, it 
might not be possible to solve analytically the algebraic equations of motion for the x's, for 
example, if the equations of motion include high-order polynomials. (It might also occur 
that the solution of the equations of motion for the x' s is n °t unique. In that case, the 
action and the system it describes splits into separate self-consistent Hamiltonian theories.) 
In Appendix A the familiar Hamiltonian form of the action for Einstein gravity coupled to 
Maxwell electrodynamics is derived using the algorithm developed in Sec. 3. 

A fourth path integral method considered by Iyer and Wald is (iv) the calculation of 
exp(S BH ) as the enhancement factor for the rate of black hole pair creation relative to the 
pair creation rate for matter distributions. They show || that method (iv) yields the entropy 
expression (1.1) when applied to an arbitrary diffeomorphism invariant theory. In Ref. [15 



it was shown that the enhancement in the black hole creation rate, method (iv), must 
agree with the entropy as calculated from the microcanonical functional integral, method 



i). Therefore the results of Ref. Jl5| along with those obtained here constitute a derivation 



of Eq. (1.2) as the enhancement factor for black hole pair creation in a diffeomorphism 
invariant theory. 

The analysis presented here applies to any stationary spacetime with bifurcate Killing 
horizon. This includes not only certain black hole spacetimes but also, for example, Rindler 
spacetime. With periodic identifications in the extra dimensions, the bifurcation surface 
7i of the Rindler horizon has the topology of a {D — 2)— torus. Equation (1.2) gives the 
associated entropy. For definiteness, I will typically use the terminology appropriate for 
black hole spacetimes. Some key results concerning the surface gravity of a bifurcate Killing 
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horizon, which are used in the analysis of Sec. 4, are derived in Appendix B. 



II. MICROCANONICAL FUNCTIONAL INTEGRAL 

In the microcanonical functional integral formalism, the density of states v is expressed 
directly as a path integral PJID|: 



v = Y.\ DgD^exp(S[g^}) . (2.1) 

M J 

Here, S is the action, which is a functional of the metric g a b and a collection of matter fields 
denoted by ip. Also, YjM denotes a sum over manifolds M of different topologies, subject 
to the requirement that the boundary dM should have topology B x S 1 . For the purpose 
of describing the thermodynamics associated with a horizon, it is most convenient to choose 
B to have the same topology as the bifurcation surface. Thus, for the case of black hole 
spacetimes, B is a (D — 2)-sphere. For the case of Rindler spacetime B is a (D — 2)-torus. 

The boundary conditions on the metric and matter fields in the path integral for the 
density of states v involve fixation of those quantities on dM that characterize the states 
of the system. In the terminology of traditional thermodynamics, these are the extensive 
variables including, for example, internal energy and electric charge. These quantities appear 
at the classical level as functions of the canonical variables q a and p a . Thus, consider the 
action S written in Hamiltonian form, 

S[X,q,p] = i dt d d x[p a q a — X A C A (q,p)^j + (boundary terms) , (2.2) 

where d — D — 1 is the dimensionality of space E. The Lagrange multipliers are denoted 
by X A and C A (q,p) are the constraints. The particular form (2.2) for the action follows from 
spacetime diffeomorphism invariance and the assumption that under reparametrizations in t 
the canonical variables transform as scalars and the Lagrange multipliers transform as scalar 
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There are two types of boundary terms that appear in the Hamiltonian form of the 
action for a diffeomorphism invariant theory on a manifold M. The first is a term on the 
boundary dM = B x S 1 of the spacetime manifold. By the argument given in Ref. [15|, the 



action (2.2) appropriate for the density of states v contains no such boundary terms at dM. 
Otherwise, if boundary terms at dM were present, the boundary conditions would include 
fixation of quantities that depend on the Lagrange multipliers (and, hence, do not depend 
solely on the canonical variables q a and p a ). 

The second type of boundary term that appears in the Hamiltonian form of the action 
arises only if the boundary of space E includes an element 7i in addition to the generic 
leaf B of the foliation of dM = B x S 1 ; that is, if <9E = TC U B. This situation occurs in 
particular when the spacetime manifold has topology M = B x M 2 and the leaves of the 
foliation terminate at a common surface Ti., considered to be the "origin" of the IR 2 plane. 
Such boundary terms are derived as follows. Start with the action in Lagrangian form, 
expressed as an integral over M (plus possible boundary terms). Now excise a region from 
M surrounding H, so that M has the product topology B x (annulus) = E x S 1 (where 
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E = B X I, with I a real line interval). The boundary DM. then consists of two copies 
of B x S 1 , where one copy coincides with the original boundary of Ai and the other copy 
coincides with the boundary of the excised region. The passage from the Lagrangian form of 
the action to the Hamiltonian form of the action proceeds as usual, with various boundary 
terms appearing at the boundary of the excised region. One then takes the limit in which 
the excised region shrinks to zero, being careful to insure that the geometry is smooth at 
7i. The second type of boundary term is a term on the boundary of the excised region that 
survives this limit. 

The entropy of a stationary black hole is computed as follows PJIDfl- First, express 
the Lorentzian solution in stationary coordinates, ds 2 = g a bdx a dx b , ip — i/j, where g a b and 
ip are t independent. Next, choose boundary conditions for the path integral (2.1) that 
coincide with the boundary values (as constructed from the canonical data of a t — const 
slice with boundary element B) of the black hole spacetime. The path integral for v will 
have an extremum in the topological sector Ai = B x IR 2 that consists of the complex 
black hole solution ds 2 = g a b dx a dx b , if) = ip. The complex black hole is obtained from the 
Lorentzian black hole by the substitution t — > —it. The t = const slices of the Lorentzian and 



complex black hole solutions coincide in the sense that their canonical data agree fli7|P,|15 
In particular, the data on the boundary element 7i of the spatial slices of the complex black 
hole coincide with the data on the bifurcation surface of the Lorentzian black hole. In the 
zero-loop approximation the density of states is given by v « exp(S[g, ?/>]), where S[g,ijj} is 
the action evaluated at the complex black hole solution. The entropy of the black hole is 
then 

S BH *S\g,$], (2.3) 

the logarithm of the density of states. 

When the action is expressed in Hamiltonian form, the evaluation of the entropy in 
Eq. (2.3) is simple. Because the complex spacetime g, ip is a stationary solution of the 
classical equations of motion, both the p a q a terms and the constraint terms in Eq. (2.2) 
vanish. The only contribution to the entropy comes from the boundary terms. As discussed 
above (see Ref. there are no boundary terms at dM.. There are, however, boundary 

terms at 7i. The entropy arises entirely from the evaluation of these terms at the complex 
black hole solution. 



III. ACTION FOR DIFFEOMORPHISM INVARIANT THEORIES 

The action for an arbitrary diffeomorphism invariant theory of the metric g a b and tensor 
matter fields ip can be expressed in the manifestly covariant form || 

S[g,ijj] — i d D xy/^gC , (3.1a) 

J M 

c = c( gab, TZ-bcdei V ' ai T^-bcdei ■ ■ ■ i • • ^ a m )T^-bcdei ^ i V Ql ^, . . . , V( 0l ' • • V a< )^] , (3.1b) 
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where V a is the spacetime covariant derivative.^ I will use the action (3.1) as the starting 
point. Note, however, that the derivation of the entropy (1.2) in this paper is not restricted 
just to the case of tensor matter fields. For example, ip can include the components A% of 
the Yang-Mills connection (a is the internal index), where the covariant derivative V a of 
Eq. (3.1b) acts on A" as a collection of covariant vectors. Likewise, ip can include the tetrad 
field (e M ) a . In this case, the Lagrangian should include a term A ab [g a b — (e fl ) a ri flu (e lJ )b] that 
links the tetrad to the metric, where A ab is an independently varied field (included among 
the ip) and rj^ = diag(— 1, +1, . . . , +1). With the tetrad appearing as a dynamical variable 
one can also include coupling to the Dirac field (see, for example, Ref. ]B|). 

The goal of this section is to place the action (3.1) in Hamiltonian form (2.2), and 
thereby derive the relevant boundary terms at TC. I will assume that the manifold topology 
is M. = E x I. In Sec. 4, where the density of states (2.1) is evaluated, the factor / is 
periodically identified to form a circle S 1 . 



A. Elimination of derivatives of the Riemann tensor 



The first step in the derivation of the Hamiltonian form of the action (3.1) is the elimina- 
tion of derivatives of the Riemann tensor. The highest derivative, namely the m th derivative, 
can be eliminated as follows. Introduce a set of auxiliary fields jj^"' ambcde and VT*... a 6c(fe , 
and write the action as 



s[g, ij, u m , v m ] = i J m d D xV=g{c(g, n, v ai n, . . . , v (ai - • ■ v am _ 1} n, v c 



ai---a m ' 



+ U, 



V ...V V — v m 



(3.2) 



Note that the indices on the Riemann tensor, and the corresponding indices on U m and 
V m , have been suppressed. Also, the notation ^'s is used for the matter fields and their 
derivatives. The actions (3.1) and (3.2) are equivalent. This is demonstrated by substituting 
the solution of the classical equations of motion for U m and V m , namely, 








5S 



1 5S 

v ai---a m 



V ... V 7? — V m 

V(ai v a m )'V- v ai --a m 

DC 



Tja\—a n 



+ 



dV m 

ai---a„ 



(3.3a) 
(3.3b) 



into the action (3.2). The result is the action (3.1). Now integrate by parts in Eq. (3.2) to 
remove one derivative from V (ai ---V am) 7?., and discard the boundary term. This leads to the 
action 



2 The work of Anderson and Torre (T^j implies that the Lagrangian (3.1b) can be written in 
terms of covariant derivatives of the Riemann tensor in which the symmetrization over covariant 
derivatives is extended to the second and fourth slots of the Riemann tensor itself. The arguments 
of this section remain valid whether or not the Lagrangian is written in this way. 
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S[g,ip,U m , V 

d D x 



J m d D x^{c(g, K, V ai K, . . . , V (ai - V^TZ, V™.. am ^s) 



(3.4) 



in which the highest derivative of the Riemann tensor is the (m — l) th derivative. 

The action (3.4) yields the same equations of motion as the original action (3.1), but (3.4) 
is not entirely equivalent to (3.1) because boundary terms were discarded in its derivation. 
The change of boundary terms implies a change of boundary conditions for the variational 
problem. However, the boundary terms at dAi that should be present in the final Hamil- 
tonian form (2.2) of the action are known: there should be no boundary terms at dAi. 
Thus, we are free to discard any boundary terms that arise through integration by parts 
in spacetime. At the end of the analysis, any remaining boundary terms at dAi must be 
eliminated from the Hamiltonian form of the action anyway. 

The algorithm described above can be iterated until all derivatives of the Riemann tensor 
have been eliminated. In the process, sets of auxiliary variables £7 m _i, F m_1 , . . ., Lq, V 1 
are introduced, which serve to eliminate the (m — l) th , . . ., (l) th derivatives of 1Z bcde . The 
resulting action is 



S[g^,U m ,V m ,...,U 1 ,V 1 \ 

= i j M d D x^gl [ C(g, n,v\...,v m , V>'s) 



(V ai U?)K + U?V a \ (v a2 ur 2 )v a \ + ur 2 v c 



r2 



(V„ U ai '" am )V 

\ v a-m^ m 1 c 



m—1 , rrsi-a m ym 

a l'" a m— 1 m Ol"-aj7 



(3.5) 



Now isolate the Riemann tensor by introducing one more set of auxiliary fields, Uq and V . 
(This is one more iteration of the algorithm, but without the integration by parts). The 
action becomes 

S[g^,U m ,V m ,...,U ,V } 

= if d D x v ^\c(g,V°,...,V m ,^s) 



+ 



U K - UnV 



0' 



(v ai unv° + urv a \ 



(X7 Tj ai "' am )V 

\Va m u m ) y c 



m—1 , Tja\---a m -\rm 

ai-a m -i ' m ai---a„ 



(3.6) 



Observe that when the V equations of motion hold, 

dC 



= -Un 



v ai c/r + 



dV° ' 



U — ~ U m-l — V a m U m 

dC 



n _ Tja\---a m 



+ 



gym 



dC 



gym-i 
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and the U equations of motion hold, 

= K - V° 



o = v ai v°-v a \ , 

n = v v™^ 1 — v m 

w v (dm a\---a m -\) r ai---a v 



the variable Uq equals 

dC _ / dC 



jjbcde " Y7 

~~ a/T5 \ V ai 



d(n bcde ) ai \d(v ai n 



-bcde i 



- +( - irv -- v -( iro^) )' (37) 

This is the variational derivative M of the Lagrangian C with respect to the Riemann tensor 

B. Elimination of higher order derivatives of the matter fields 

The second step in the derivation of the Hamiltonian form of the action is the elimination 
of all but the first covariant derivatives of the matter fields. This can be achieved by 
applying the same algorithm that was used in the elimination of derivatives of the Riemann 
tensor. The resulting action depends on the first covariant derivative of the following fields: 
Ui, . . . , U m , some (or all) of the matter fields if), and some of the auxiliary fields that were 
introduced in the elimination of the higher-order derivatives of ip. I will denote these 
fields collectively by 1 J r '. The covariant derivatives V a \l/' can be isolated as linear terms in 
the Lagrangian through the introduction of yet another set of auxiliary fields, just as the 
Riemann tensor was isolated by the introduction of the fields Uq, V°. The action now takes 
the form 

S[g, V] = i[ d D x^{ut d n abcd + f(g, V *')} , (3.8) 

where V a $' appear linearly in / with coefficients that are independent variables. In Eq. (3.8), 
^ denotes the original matter fields ip, the auxiliary fields U and V, and the auxiliary fields 



3 It is possible to carry out the analysis above without the fields V. For example, in order to 
eliminate the m th derivative of the Riemann tensor from the action (3.1), one can perform a 
Legendre transformation in which V(a m (V ai - ■ ■ V 0m _ 1 )72.) play the role of velocities and U^"' am 
play the role of momenta, and follow this with an integration by parts. (Equivalently, V m can be 
eliminated from the action (3.4) by substitution of the solution of the V m equation of motion.) 
One must allow for the possibility that the relationship between the velocities and momenta is not 
invertible, signaling the presence of constraints. The key equation (3.7) and the form (3.8) of the 
action (see below) can be deduced in this way, in spite of the fact that the constraints are not 
known explicitly. 
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that were introduced in the elimination of the higher-order derivatives of ip and the isolation 
of the first derivatives V a \l/'. Thus, ty' is the subset of fields \1/ that appear differentiated 
in the action. The presentation below is simplified if we assume that the fields are 
covariant in their tensor indices. (As discussed at the beginning of this section, some matter 
fields might also carry internal indices, such as a Yang-Mills index or a tetrad index.) This 
assumption entails no loss of generality — for each field with a contravariant tensor index 
that appears differentiated in the action, say ip a , a simple change of variables if) a = g a btp b , 
9ab = 9ab allows us to replace ip a with ip a as the fundamental variable. 

C. Space time decomposition 

The third step in the derivation of the Hamiltonian form of the action is the introduction 
of a space-time split. Let spacetime have topology Ai = H x I and let t label the hyper- 
surfaces of the foliation E. The unit normal of the hypersurfaces is u a = — iVV a t, where 
N = [— (V a t)(7 a6 (V;,t)] _1 / 2 defines the lapse function. The hypersurface metric is defined by 
h a b = Qab + u a Ub, so the spacetime metric g ab becomes 

g a b = Kb - u a u b . (3.9) 

The Gauss, Codazzi, and Ricci equations imply (see, for example, Ref. p0| ) 

Tlabcd = Rabcd + 2K a[c K d]b + 4(D [a K b][c )u d] + A(D [c K d][a )u b] 

-4u la (£ u K b]lc + K e b] K e[c + (D b] D [c N)/N)u d] , (3.10) 

where £ u is the Lie derivative along u a and R a bcd, K ab = — \£ u h ab , and D a denote the 
Riemann tensor, extrinsic curvature, and covariant derivative of the t = const hypersurfaces, 
respectively. 

Under the space-time decomposition the tensor fields \1/ are projected normally and 
tangentially to the t = const hypersurfaces. For example, ip a is split into its projections 
u a i\) a and h a ip b and tp a is split into u a tp a and h^t\} h . The first derivative V a $' of a covariant 
tensor is split into hypersurface covariant derivatives and normal Lie derivatives of the 
projections of ^' . For example, for a scalar field ip we have V a ip = —u a £ u ip + D a %p and for 
a covariant vector field ip a we have 

V a ^ b = u a u b £ u {u c ^) c ) - u a £ u {h c b ip c ) + D a {h c b ip c ) - u b D a (u c ip c ) 

+K ab u c ip c - 2u (a K c b) ij c - u a u b ijj c (D c N)/N + u a u c ^ c {D b N)/N . (3.11) 

Note that the Lie derivatives £ u and the spatial covariant derivative of the lapse function, 
in the combination (D a N)/N, appear linearly in V a ip b . Also note that £ u (h b ip c ) is a spatial 
tensor; that is, u b £ u (hlip c ) = 0. The decomposition for the derivatives of higher rank 
covariant tensors is similar to that in Eq. (3.11). I will use V^f as a shorthand notation for 
the normal and tangential projections of the fields 

With the space-time split described above, the action (3.8) becomes 
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\j M d ' \ A 1 



R abcd + 2K ac K bd + 8U£ bcd u d (D a K bc ) 



-AU« bcd u a u d [£ u K bc + K e b K ec + (D b D c N)/N 
+f(h ab , h ab , VV, £ V (V*% D a (V*'), K ab , (D a N)/N)} , (3.12) 

where £ u (V^f) and (D a N)/N appear linearly in /. The action (3.12) contains a second 
time derivative in the term £ u K bc = —\£ u £ u h ab . This can be removed by promoting the 
extrinsic curvature K bc to an independent variable. Thus, introduce an auxiliary variable 
P ab and write the action as 



NP 



ab 



£ u h a b + 2K ab 



+ Ut d [Rabcd + 2K ac K bd ] + 8^ Ut d u d (D a K bc ) 
- A^)U* bcd u a u d [£ u K bc + K e b K ec + (D b D c N)/N 

+ V=9 f(h ab , h ab , VV, £ U {VV), D a (W% K ab , (D a N)/N) 



(3.13) 



The action (3.12) is recovered when the solution of the P ab , K ab equations of motion is 
substituted into the action (3.13). 

Now choose a time flow vector field t a such that t a V a t = 1, and define the shift vector by 
V a = h b l t b (not to be confused with the auxiliary fields V°, . . ., V m ). The Lie derivative £ u , 
acting on h ab , K ab , and the covariant tensors V^', is expressed as N£ u = £ t — £y. With 
the fields mapped from Ai to E x J, the Lie derivatives £t along t a become ordinary time 
derivatives (denoted by a dot) and ^f^g = Ny/h where h is the determinant of the metric 
hij on E. The action (3.13) becomes 



S[N,V,h,K,P,VV] 



=i I dt L 



d d xl P ij 



h i:j - + 2NKij 



+NVhUi jke [R ijke + 2K ik K j£ ] - 8NVhUt u (D t K jk ) 
-4VhU™ ju [k i:j - £ v K ij + NKj°K kj + (DiDjN) 

+NVhf(hii, h ij , VV, £ U (W), Dip*'), Kij, (DiN)/N) \ , 



(3.14) 



where N £ U {V S S>') = (V^')' — £v(P^') and i, j, ... are indices for tensors on E. Also, the 
notation C/ abcu = -U$ bcd u d and U$ bcu = U$ bcd u a u d has been used. 



D. Hamiltonian form of the action 

The action (3.14) is linear in time derivatives and Lagrange multipliers. That is, each 
term in the Lagrangian of Eq. (3.14) depends linearly on either the time derivative of a field, 
the lapse function N, or the shift vector V 1 . (Recall that £ U (V^') and (D a N)/N appear 
linearly in /.) Therefore each term in the Lagrangian transforms as a scalar density under 
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reparametrizations in t, and no terms such as V l hijV^ /N appear. Moreover, the coefficients 
of hij, Kij, and {V^')' are independent variables. Specifically, the coefficient of is P*- 7 , 
the coefficient of is —A\fhUQ %:)U , and the coefficient of (V^ 1 )' is given in terms of the 
auxiliary fields that were introduced in the process of isolating V a ^' as a linear factor in the 
action. It follows that these coefficients (denoted p a ) are the canonical momenta conjugate 
to the coordinates h^, K^, and V§>' (denoted q a ). After all, consider what happens if one 
tries to identify both q a and p a as coordinates, and to define conjugate momenta n^, and IT?. 
The definition of these momentum variables leads to sets of second class constraints, namely, 
H"a = p a and lip = 0. Elimination of these constraints through the Dirac bracket effectively 
eliminates the new momenta, and reveals the interpretation of q a and p a as canonically 
conjugate variables. 

Now remove the spatial derivatives from the lapse function N and the shift vector V 1 
through integrations by parts. The action (3.14) takes the "almost Hamiltonian" form 

S[N, V, q,p, X ]=tJ dt jf d d x(p a q a - NC ± (q,p, X ) - V%{q,p)) 

+i J dt J ^xy/a^-AmU^DjN + (terms ~ N and V*)) , (3.15) 

where n l is the outward pointing unit normal of the boundary of space, <9£, and a is the 
determinant of the metric on <9£. The "terms ~ N and V in that appear at the boundary <9£ 
are proportional to the undifferentiated lapse function N and undifferentiated shift vector 
V 1 . Observe that the action for any diffeomorphism invariant theory can be put into the 
form (3.15), since no special assumptions were made in its derivation. 

The action (3.15) is not quite in Hamiltonian form because it contains certain extra 
undifferentiated variables X i n addition to the canonical variables q a and p a and Lagrange 
multipliers N and V 1 . These variables include, for example, the normal and tangential 
projections of the auxiliary fields V°, . . ., V m . Note that the x's appear only in the Hamil- 
tonian constraint C±. They do not appear in the momentum constraint Cj, or in the boundary 
terms proportional to V\ because Cj and the boundary terms proportional to V 1 all orig- 
inate from the Lie derivatives £ u in the combination p a (q a — £yq a )- The x' s also do not 
appear in the boundary terms proportional to N. To see this, one should recall that these 
boundary terms arise through integration by parts that eliminate spatial derivatives of N. 
Inspection of Eq. (3.10) shows that the spatial derivatives of iV appear in the combination 
N£ u K bc + D b D c N. Thus, the coefficient of DiDjN in the action equals the momentum 
conjugate to K^, and involves no x' s - Likewise, inspection of Eq. (3.11) shows that the 
spatial derivatives of iV appear in combination with Lie derivatives in such a way that the 
coefficients of DiN involve the momenta conjugate to V^' and V^' itself, but no x's- 

A true Hamiltonian form of the action is obtained from Eq. (3.15) by elimination of the 
variables x through the solution of their algebraic equations of motion, 

% 5S dC± 

^^ = ^ = °- 3 - 16 

N5 X d X K ' 

That is, one solves the set of equations (3.16) and inserts the solution back into the action. 
If the matrix of second derivatives of C± with respect to the x's has vanishing determinant, 
then the equations (3.16) are not independent and cannot be solved for all of the x' s as 
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functions of the canonical variables q a and p a . In that case the equations (3.16) include 
constraints on the canonical variables. In any case, when the conditions (3.16) are imposed 
the Hamiltonian constraint Cj_ is independent of x- 

The constraints that arise through the elimination of the x's must be incorporated into 
the action principle via Lagrange multipliers. Let C A denote the complete set of constraints 
for the system — the constraints that arise through equations (3.16) as well as the Hamilto- 
nian constraint C± and momentum constraints C{. Likewise, let X A denote the complete set 
of Lagrange multipliers associated with the constraints C A , including the lapse function N 
and shift vector V 1 . The action, which is now in Hamiltonian form, reads 

S[X,q,p]=i J dt Jj d x(p a q a - \ A C A {q,p)) 

+i J dtj ^xy/a^-AmU^DjN + (terms ~ N and V*)) . (3.17) 

This is Eq. (2.2) with the boundary terms displayed somewhat more explicitly. 

The details of the elimination of the variables x must be carried out on a case-by-case 
basis. The instructive example of Einstein gravity coupled to Maxwell electrodynamics is 
presented in Appendix A. 



E. Comments 



Several comments are in order. First, consider the situation in which one of the con- 
straints, say, Ci, is simply one of the momentum variables, say, p\. Then the equation of 
motion for the Lagrange multiplier Ai is p\ — 0, and the equation of motion for p\ yields an 
expression for Ai in terms of the other Lagrange multipliers, q a , and p a . These equations 
can be used to eliminate p\ and A 1 from the variational principle — effectively one just sets p\ 
equal to zero. In this way the pair q 1 , p 1 is removed from the list of dynamical variables in 
the theory, although in general the action still depends on q l (undifferentiated in time). The 
resulting situation is similar to that encountered in the "almost Hamiltonian" form (3.15) 
of the action, in that the action depends on an extra variable. (The key difference is that 
the constraints might depend on spatial derivatives of q 1 , and also that q l might appear in 
the boundary terms. The variables Xi on the other hand, appeared only undifferentiated in 
the Hamiltonian constraint.) Now one can attempt to eliminate q l through the solution of 
its equation of motion. If the q 1 equation of motion can be solved for q 1 , then insertion of 
this solution into the action yields a new Hamiltonian form of the action for the system in 
which q 1 and p\ are completely excluded. If the q 1 equation of motion depends only on the 
canonical variables (other than q 1 and pi), then q 1 is a Lagrange multiplier and should be 
left alone. It might happen that the q 1 equation of motion cannot be solved for q 1 , but also 
does not yield a constraint. In this case one can always stick to the Hamiltonian form of the 
action that includes q 1 , p\ and the constraint C\ —p\. 

Although the situation in which one of the constraints is equal to a momentum variable 
might appear to be of academic interest only, it in fact occurs in the examples of Maxwell 
electrodynamics and Einstein gravity. In electrodynamics the variable that plays the role 
of q 1 is the normal projection of the electromagnetic potential, which becomes a Lagrange 
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multiplier for the Gauss's law constraint. In Einstein gravity one must first perform a 
canonical transformation on the variables q a , p a to bring the action to a form in which one 
of the momentum variables is constrained to vanish. The variables (q 1 and pi) that are 
eliminated in this way are the extrinsic curvature and its conjugate. Details can be 
found in Appendix A. 

As a final comment, observe that the Hamiltonian H — J E d d x\ A C A + (boundary terms) 
obtained from the action functional (3.17) is not necessarily either the total Hamiltonian or 
the extended Hamiltonian [[Tj|. If the extended Hamiltonian of the system is desired, one can 
start with the Hamiltonian H and treat the constraints C A = as primary constraints. The 
preservation in time of the primary constraints can lead to secondary constraints. One then 
proceeds to the classification of constraints as first or second class, and to the construction 
of the extended Hamiltonian |]I6 |. For the purpose of this paper, it is not necessary that 
the extended Hamiltonian appear in the Hamiltonian form (3.17) of the action. What really 
matters is that the variational principles based on the action functionals (3.1) and (3.17) 
are equivalent. 



IV. BLACK HOLE ENTROPY 

We are now in a position to compute the entropy of a stationary spacetime with bifurcate 
Killing horizon using the microcanonical functional integral method. It is assumed that the 
spacetime metric and matter fields {g a b,ip} satisfy the classical equations of motion that 
follow from the action (3.1). Let t a denote the Killing vector field that vanishes on the 
bifurcation surface 7i. Consider a spacelike hypersurface S whose boundaries consist of an 
outer boundary B and the bifurcation surface 7i. Note that S lies within a single "wedge" of 
the spacetime where t a is timelike. Now extend Eo into a foliation of the wedge by stationary 
hypersurfaces t = const, where t a V a t = 1. Also choose t a as the time flow vector fieldfj 
Then the solution {g a b, i>} can be written in Hamiltonian form {A, q, p}, where A, q, and p 
are t-independent. 

As discussed in Sec. 2, the entropy is obtained by evaluating the action (2.2) at the 
complex solution {A, q, p} with periodic identification in t. The complex solution {X,q, p} 
is obtained from the real Lorentzian solution {X,q, p] by the substitution t —>■ —it. Under 



4 One can choose a time flow vector field t a — £l<j) a where fi is constant and 4> a is a spatial Killing 
vector field, if such a Killing vector field exists. This changes the shift vector by — £l(ft a , so the 
argument given in Appendix B that V a vanishes at TC no longer holds. However, the overall results 
of the analysis are unchanged because in the action (4.1) below the extra nonvanishing boundary 
terms at H just cancel corresponding boundary terms at B. This can be seen by reversing the steps 
that generate the boundary terms proportional to V a : Sum the identity p a £yq a = over canonical 
pairs, integrate over S, then integrate by parts and use the momentum constraint. Experience with 
black hole thermodynamics in the context of Einstein gravity [1^,21] shows that the shift vector 
defined as the spatial projection of t a , not t a — £l(j) a , is the physically correct definition for the 
product of inverse temperature and chemical potential. 
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reparametrizations in t the canonical variables transform as scalars and the Lagrange mul- 



tipliers transform as scalar densities, so it follows |15[ that A = —iX, q = q, and p = p. The 
orbits of the Killing vector field t a in the complex spacetime form closed curves (circles) 
around the bifurcation surface TC. 

According to the discussion of Sec. 2, the action (2.2) contains no boundary terms at 
the spacetime boundary dM. = B x S 1 where the boundary data appropriate for the micro- 
canonical functional integral are fixed [[TJJ. The boundary terms at TC are just the boundary 
terms displayed in Eq. (3.17). Thus, the action takes the form 

S[X,q,p]=iJ dt Jj d x(p a q a -X A C A (q,p)) 

+i f dt f d^xVaf-AmU^DjN + (terms ~ N and V*)) , (4.1) 

where <9X = TC U B. The boundary terms at TC should be understood in terms of a limiting 
procedure, as discussed in Sec. 2. That is, the boundary terms are defined by an integral 
over the boundary of an excised region that surrounds the bifurcation surface TC, and the 
limit is taken as the excised region shrinks to TC. 

In Appendix B it is shown that the lapse function N and shift vector V a , and hence also 
N and V a , vanish in the limit as the bifurcation surface is approached. Thus, all boundary 
"terms ~ N and V in in the action (4.1) vanish upon evaluation at the complex solution. 
This assumes that the coefficients of the "terms ~ N and V 1 " , which depend solely on the 
canonical variables, are well behaved at the bifurcation surface TC. Now, because {A, q, p} is 
a stationary solution of the classical equations of motion, the p a q a terms and the constraint 
terms vanish in the evaluation of the action. The result is that the entropy (2.3) becomes 

S BH « S[X,q,p] = -Ai [ dt [ d^x^riiUo^DjN , (4.2) 

The right-hand side of this expression is evaluated at the complex solution {X,q, p}. How- 
ever, for notational simplicity, the bars have been omitted from the p's and q's in Eq. (4.2). 
No ambiguity arises since the canonical variables for the Lorentzian and complex solutions 
agree. The bar is retained on the lapse function since the Lagrange multipliers for the 
Lorentzian and complex solutions differ by a factor of —i. In the analysis below I will 
continue the practice of placing bars or tildes only over the Lagrange multipliers. 

As shown in Appendix B, the gradient of the lapse function is related to the surface 
gravity R by YmvDiN = — limkni, where the limit is taken in which the bifurcation surface 
TC is approached along a t = const hypersurface S. (Note that n l is the outward pointing 
normal of <9£ at the bifurcation surface, so n l points "radially inward" towards TC.) Thus, 
we obtain 

S B h ~ 4 / dt [ d^xy/amU^nik . (4.3) 
Js 1 Jh 

Since the surface gravity of a spacetime with bifurcate Killing horizon is constant over the 
horizon Jffi] , k can be removed from the integral over TC. Now, the proper circumference 
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of the circular orbits of t a is J s i dt\J —N 2 + V^V^.Q From Appendix B we have the result 
limV/iV = 0, so the proper circumference, in the limit as the bifurcation surface is ap- 
proached, equals J s i dtN . The expression k = —limn^DjN then shows that J s i dtk equals 
the rate of change of circumference with respect to radius for these orbits. The complex 
geometry will be smooth at 7i, and satisfy the classical equations of motion there, only if 
the period in S 1 is chosen such that J s i dtk = 2tt. The entropy is then 

S BH « 8tt f ^xy/dmU^rij . (4.4) 

The right-hand side of this expression for S BH depends only on the canonical variables, so 
it can be evaluated either at the complex solution {A, q, p} or at the Lorentzian solution 

Recalling the definition [/q 6cu = U^ cd u a Ud and using the expression l a b = 2 lim u {a n^ for 
the binormal of 7i, we have 

S BH « -2tt / tf-^xyfdlaljjt* . (4.5) 

This is the main result, Eq. (1.2), for the entropy of a spacetime with bifurcate Killing 
horizon. Here, UQ bcd is the variational derivative (3.7) of the Lagrangian with respect to the 
Riemann tensor. 



V. OTHER PATH INTEGRAL METHODS 
A. Hilbert action surface term 



The relationship between the Hilbert action surface term method (ii) [|iT| , |12| and the 
microcanonical functional integral method (i) can be understood as follows. Consider first 
the logical outline of the microcanonical functional integral method. The action S is the 
integral of the Lagrangian £ over the manifold Ai = B x IR 2 . The integral is split into two 
pieces, an integral over the excised region B x disk that contains the bifurcation surface 7i, 
and an integral over the remainder of the manifold, B x annulus. Schematically, we have 



S= / C 

Jm 2 



I C+fc, (5.1) 
Jd J a 



where D and A refer to the disk and the annulus, respectively, and the factor B has been 
suppressed for notational simplicity. The first integral vanishes in the limit in which the 
excised region D shrinks to 7i, under the assumption that £ is smooth. The second integral 
is written in Hamiltonian form, which yields a volume integral S c [the integral of p a q a — \ A C A 
in Eq. (4.1)] and terms at the "inner" boundary dAi and "outer" boundary dA Q of the 



5 This is / ds, where ds 2 = -N 2 dt 2 + h ij (dx i + V i dt)(dx j + V j dt) with dx i = 0, evaluated at the 
complex solution. 
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annulus. The terms at the outer boundary dA Q are discarded. In the limit as the disk 
shrinks to Ti,, the terms at the inner boundary dAi become boundary terms BT^ at the 
bifurcation surface [the boundary terms in Eq. (4.1)]. Thus, the action becomes 

S -> S c + BT H . (5.2) 

When S is evaluated at the complex solution, the canonical action S c vanishes so that 

<S — > BT n . (5.3) 

Only the term proportional to the gradient of the lapse function N remains in BT-^ since, 
as shown in Appendix B, iV and V 1 both vanish at the bifurcation surface. 

In the Hilbert action surface term method (ii), the action integral is split according to 

S= ( C= f C+ f £ + I £+ I i, (5.4) 

JIR 2 JD JdD J A JdAi 

where J I is the "Hilbert action surface term"; that is, j £ is the surface term that must 
be added to / £ such that the boundary conditions include fixation of the metric on the 
boundary. Equation (5.4) is, of course, equivalent to Eq. (5.1) since the Hilbert action 
surface terms at dD and dAi cancel one another. As in the microcanonical functional integral 
method, the integral over D vanishes in the limit in which the excised region shrinks to Ti. 
Also the integral over A can be written in Hamiltonian form. The resulting action is 

£ + S c + BT' n , (5.5) 

dD 

where BT^ = BT^ + j aA . £. (Again, the terms at the outer boundary are discarded.) 

The integrand in BT' U must be linear in the lapse and shift in order to transform prop- 
erly under reparametrizations in t. But, in fact, the integrand in BT' n cannot depend on 
spatial derivatives of N or V 1 . (Actually, spatial derivatives in a direction tangent to the 
boundary are allowed, since these can be removed through integration by parts.) This can 
be understood as follows. The boundary conditions appropriate for S c + BT^ (which equals 
J A £ + j dA . £ plus terms at the outer boundary dA ) include fixation of the induced metric 
on dAi, by definition of the Hilbert action surface term. Therefore, with the induced metric 
on dAi denoted by 7 mn , we have 

5{S C + BT' n ) = (eom's) + (bt's at dA Q ) + / ii mn 5 lrnn + (other bt's at dA { ) (5.6) 

for some n mn . Here, "eom's" are terms that yield the classical equations of motion and 
"bt's at dA Q n are boundary terms at dA Q . The "other bt's at dA" are boundary terms at 
dAi that involve variations of various matter fields and auxiliary fields, but do not involve 
variations of 7 mri or variations of derivatives of / ~f mn - The induced metric 7 mn on the (D — 1 
dimensional) surface dAi can be split into a lapse, shift, and (D — 2 dimensional) spatial 
metric using the slices t = const and time flow vector field t a induced on dAi. If the slices 
t = const are orthogonal to dA { , so that the unit normal of the slices lies in dA i: then 
the lapse and shift components of 7 m „ are just the restrictions of N and V % to dAi. If the 
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slices t = const are not orthogonal to dAi, then the lapse and shift components of 7 mn are 



constructed algebraically from N and V 1 through simple kinematical boost relations |f23 
Consequently, the boundary terms at dAi in Eq. (5.6) depend on the variations of N and 
V 1 , but not on the variations of their derivatives. Since S c contains no derivatives of iV or 
V 1 , Eq. (5.6) shows that BT' n cannot contain derivatives of N or V 1 . 

When the action of Eq. (5.5) is evaluated at the complex solution, the canonical action S c 
vanishes. The boundary term BT' n is zero since its integrand is linear in the undifferentiated 
lapse and shift, and the lapse and shift vanish at 7i. Therefore 



S - / £, (5.7) 

JdD 

which shows that the entropy S BH ~ S(X,q,p) equals the Hilbert action surface term for 
a small disk surrounding the bifurcation surface 7i, evaluated at the complex solution. In 
effect, what has been shown is that the Hilbert action surface term must include the negative 
of the particular boundary term displayed in Eq. (4.1) that is proportional to the gradient 
of the lapse function. The minus sign is compensated by the fact that the normal n l of dD 
points away from 7i. 



B. Conical deficit angle 

The starting point for the canonical deficit angle method (iii) |T^JT3| j6]] for computing 
black hole entropy is the (grand canonical) partition function 

Z(P) = W DgDi, expfebj,^) ■ (5.8) 

M J 

Here, Sp is the Hilbert action Sp = J M C + J dM £ and the inverse temperature is defined 
by the lapse component of the boundary metric 7 mn according to /3 = i J (lapse). When 
evaluated at the complex black hole solution, f3 is the proper length in the S 1 direction as 
measured orthogonally to the stationary time slices in dM.. The relationship between Z[f3\ 
and the density of states (2.1) is spelled out in detail in Ref. ||. For our present purposes it 
is sufficient to note that the density of states is a function of the internal energy E, where 
E is defined by the variation of the Hilbert action with respect to the lapse component of 
7 mn . Thus, the microcanonical action S of Eq. (2.1) and the Hilbert action Sp differ by 
boundary terms that include a term of the form —f3E. The partition function is then given 
by the Laplace transform of the density of states:^ 

Z(J3) = J ' dEu(E)e- pE . (5.9) 



To be precise, the inverse temperature is a function on the system boundary B and, correspond- 
ingly, the energy is a surface density f[7]||. Thus, the notation used here should be viewed as 
schematic. On the other hand, (3 and E can be interpreted as the "zero mode" parts of the inverse 
temperature and energy surface density. The final result (5.15) for the entropy can be understood 
in this way. 
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With the relationship v{E) exp(S BH (E)), the leading order approximation to the integral 
in Eq. (5.9) is 

\nZ(f3)^S BH (E*)-(3E* , (5.10) 
where E* is the function of (3 such that E = E*(/3) extremizes the exponential: 

dS BH {E) 



dE 



(3. (5.11) 



Equations (5.10) and (5.11) just express \nZ(P) (which is —f3 times the free energy) as a 
Legendre transform of the entropy S BH {E). 

From the relationships above it is easy to show that the entropy is given by 

S BH (E*)n\nZ(P)-p dln g Z p {P) . (5.12) 

The zero-loop approximation to the path integral (5.8) yields lnZ((3) ss S^g,^}, where g, 
■ijj is the complex black hole solution that extremizes Sp\g,i/)] among configurations whose 
proper length (period) in the S 1 direction equals (3 at the boundary dM. Then the entropy 
can be written as 

S BH (E )« Sp[g,ip] -/3 km ■^ 7 — , (5.13) 

where the limit is taken in which f3' — > (3. It is important to recognize that in taking the 
derivative with respect to /3, one does not introduce a conical singularity in the metric g 
(contrary to claims made in the literature). Rather, when f3 is varied, the parameters of 
the solution g, ip (notably the black hole mass parameter) vary in such a way that the 
configuration remains a smooth solution of the classical equations of motion. Thus, in 
Eq. (5.13), the extremal configuration with period f3' is the smooth solution denoted g', ■?/>'. 

Equation (5.13) can be evaluated with the following trick. Since the action Sp is sta- 
tionary at g', ip', one can distort this configuration without affecting the value of the action 
to first order. Thus, replace g' with the metric g v and leave ip' alone. In principle g y could 
be any metric obtained from an infinitesimal variation of g' . Consider as a particular choice 
for g v the metric whose components are identical to the components of g in the stationary 
coordinate system, but with the period in coordinate time t adjusted so that the proper 
length in the S 1 direction at dM. is f3' rather than (3. Thus, g y is a smooth, regular solution 
of the classical equations of motion everywhere except at the bifurcation surface H. At Tt 
the regularity condition J s idtk = 2n does not hold for g v , indicating the presence of a 
conical singularity. 

With the replacement of g' by g v , the entropy becomes 

<Sbh(e*) ~ sm-^ S ^Z* M ■ ( 5 - 14 ) 

(The dependence on matter fields ifj has been dropped for notational simplicity.) Now split 
the action Sp[g] into an integral (J D C)[g] over the disk plus "other terms" that consist of an 
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integral over the annulus A and boundary integrals at dAi. The integral over D vanishes 
in the limit as the disk shrinks to TC, since £ is smooth. The other terms contain integrals 
over t and are proportional to (3. Likewise, split Spi[g^} into (f D £)[g v ] plus "other terms". 
The integral over D does not vanish in this case due to the conical singularity in g v . The 
other terms are identical to the other terms from Sp[g] with the exception that they are 
proportional to f3 f . Therefore all "other terms" in expression (5.14) cancel, and the entropy 
becomes 



S BB (E*) *t -PKm 



13' -13 \Jd 



C 



(5.15) 



The entropy is thus expressed in terms of a spacetime g w with a conical singularity. The 
limit in Eq. (5.15) is taken in which D shrinks to TC and (3' approaches f3. 

The correctness of the result (5.15) can be verified as follows. Let C be the integrand 
of Eq. (3.8) .Q The term proportional to / vanishes in the limit as the disk shrinks to TC. 
Likewise, most of the terms in UQ bcd TZ a b c d do not contribute to the entropy — only the term 
that contains the curvature in the surface orthogonal to TC survives as D shrinks to TC. 
This term captures the curvature of the conical singularity. Inserting a projection onto the 
binormal of H, we have 



D 



l - [ d D x(^U^ h e ef e gh e ab e cd TZ abcd ) 

4 J DxB \ J 



(5.16) 



The binormal can be written as e a b = 2u {a n^ where n a and u a are orthogonal to TC and to 
each other, and n a lies in a t = const surface. Note that u a = — iW a £ is imaginary when 
evaluated at a complex spacetime, so — iu a is the real unit vector with square +1 for the 
metric g v . The basic interpretation of the Riemann tensor gives 



n a (-iu b ' 



n c TZ 



abc 



5n d /A 



D 



(5.17) 



where 5n d is the change in the vector n c as it is parallel transported around the perimeter of 
the disk D (first in the n a direction, then in the — iu b direction) and Ad is the area of the disk. 
If the disk has deficit angle a, then 5n d = a(—iu d ). Thus, we find e ab e cd 7Zabcd = — ^a/A D 
and 



C 



D 



( d D x(^e ab e cd US bcd a/A D 

J DxB V 



(5.18) 



Now, the deficit angle is given by a/(2n) = {(3 — (3')/f3, so the entropy from Eq. (5.15) 
becomes 



S BH {E*) w -2tt ^d D - 2 x(^e ah t cd XJ, 



jabcd \ 
° ) 



(5.19) 



7 Iyer and Wald || conjectured that the canonical deficit angle method (iii) is limited to theories 
in which the Lagrangian is a linear function of the Riemann tensor. The results here show that 
any theory of the form (3.1) can be treated by the conical deficit angle method if the action is first 
put into the form (3.8) in which the Lagrangian is linear in 1Z a bcd- 
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In the integrand above, e ab e c dUQ bcd is a spacetime scalar and in particular it is invariant under 
reparametrizations in t. Consequently the integrand can be evaluated at the Lorentzian black 
hole solution g rather than g. The result agrees with Eq. (4.5). 
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APPENDIX A: MAXWELL ELECTRODYNAMICS AND EINSTEIN GRAVITY 

Consider Maxwell electrodynamics in D = 4 spacetime dimensions coupled to the grav- 
itational field. The electromagnetic field contribution to the action is 

S M = ~\ d A x^g-V [a A b] g ac g bd V [c A d] . (A.l) 

47T JM 

In the form of Eq. (3.8) the action becomes 

S M = i f M d 4 x^g{n ab (V a A b - A ab ) - (A lab] g ac g b % cd] )/(4n)} , (A.2) 

where V ' a A b appears linearly. The auxiliary fields are IT 16 and A ab . The "almost Hamilto- 
nian" form [cf. Eq. (3.15)] of the action is 

S M = iJ dt ^d 3 xVh{-n uu A u + U m A t - NCf - VXf) 

+i J dt (Pxy/am {NA u U ui - NA l U uu - V'Ajtt *} , (A.3) 

where the electromagnetic field contribution to the Hamiltonian and momentum constraints 
is 

Cf = —D i {W n A i ) + IF U A UU + A(n u M u ) - Yi m {KlA 3 + A m ) 
+ff u (A4 - KfAj - A m ) - U ij (piAj - KijA u - Kj) 

+(A {ij] h ik h je A m )/(47r) - (A m h ij A m )/(2ir) , (A.4a) 
Cf = -U nu DiA u - AiDjUy - 2U nj D [j A i] . (A.4b) 

The action is a functional of the coordinates (the g a 's) A u and Ai, the momenta (the p a 's) 
— \fhXl m and \^hll n \ and the extra variables (the x's) IT U , n jJI , A iu , A ui , A^-, and A uu . 

The x's appear in the action undifferentiated and only in the function Cf. In particular 
Il m , n jJ , Aj U , A U j and A^- appear quadratically in Cf and can be eliminated by the solution 
of their algebraic equations of motion. Those equations are straightforward to derive, and 
the solution is 
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rr u = -n ui , 

ITJ = -hVD [k A t] /{2Tt) , 
A-iu = DiA u — KfAj , 



A, 



— DiAj KijA u . 



Inserting this result into C™, one obtains 

Cf = -D^U^Ai) + n uu A uu + A u D t U m + 2vrir%L™ + D^D^/^Tr) 



(A.5a) 
(A.5b) 
(A.5c) 
(A.5d) 
(A.5e) 



(A.6) 



for the electromagnetic field contribution to the Hamiltonian constraint. 

The variable A uu appears in the action (A. 3), (A. 4b), (A.6) as a Lagrange multiplier 
associated with the constraint n uu = 0. As discussed in Sec. 3, the variables A uu , n uu can be 
eliminated through their equations of motion, which amounts to setting Il uu equal to zero. 
The coordinate A u conjugate to — \fhTV m remains as an extra variable in the action, which 
now reads 

S M = i J dt J^xVh^Ai - V^-AiDjU^ + 2W j D [i A j] 

-N[A n DiU ui + 27rn u %rr + D^A^A^/iAn)] J 
+i J dt J d 2 x^m {iVn u M u - V j AjU ui } , (A.7) 

The equation of motion for A u yields the Gauss's law constraint DiU m = 0. Thus, A u is 
a Lagrange multiplier. Now make the changes of variables E % = \fhli m for the momentum 
conjugate to A,- L and A t = A a t a = — NA U + V % Ai for the Lagrange multiplier. The result is 



+ A t 



S M = ij dtJ^x^Ai - V^&D^Aj 

-N[2n£ l £ t /Vh + VhD {i A 3] D [i A ]] / (An) 
-i J dt J d 2 x^fa A t n i £' 1 /Vh , 

the Hamiltonian form of the action for the electromagnetic field coupled to gravity. 
Now consider Einstein gravity in D = 4 spacetime dimensions, 

S E = i d*xyf ZI gg ac g bd "R,abcd , 

J M 



(A.8) 



(A.9) 



where Newton's constant equals l/(167r). In the form of Eq. (3.6), or equivalently (3.8), the 
action becomes 



S E = i j M d 4 xy^{ut d (n abcd - V° bcd ) + g ac g bd V» bcd ) 



(A.10) 



It is assumed that UQ bcd and V® bcd have the same symmetries as TZ a bcd- The space-time split 
leads to the action of Eq. (3.13), where the function / is given by 
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/ = h ac h bd V° bcd + 2h°»Vl bu - U^ d hlh{hlh h X fgh 

+Aus bcu h d a KHv d ° efu - ^hXv u ° cdu . (A.n) 

By mapping the fields from Ai to E x / and integrating by parts to remove spatial derivatives 
from the lapse and shift, one obtains the action in "almost Hamiltonian" form [cf. Eq. (3.15)]: 

S E = i J dtj^d 3 x{p^h t3 + Q^Kij - NCf_ ~ y'Cf) 

+i J dt J d 2 x(^/Vh)n l {-2P ij Vj - 2Q ij K jk V k + Q i] D 3 N - ND 3 Q l] ) . (A.12) 

Here, the notation Q lj = —Ay/hUQ lJU has been introduced, and the gravitational contribution 
to the Hamiltonian and momentum constraints is 

CI = -DiDjQ ij - 2P^K tJ - Q l 'K k K kj - ^fhUf l {R im + 2K ik K jt - V? m ) 

+AVhUt u mK jk - V° ku ) - VhhVvSu ~ (Q iJ + 2Vhh»)V^ u , (A.13a) 
Cf = -2D, Pi + Q jk DiK jk - 2D^ k K kl ) . (A. 13b) 

The action is a functional of the coordinates (the q a, s) hij and K ij: the momenta (the p a 's) 
Pv and and the extra variables (the x's) U l jk£ , U l jku , V° jke , V° jku , and V£ ju . 

The x variables U$ U , U^ ku , V^ ke , and V^ kn can be eliminated by the solution of their 
algebraic equations of motion. That solution is 

jjijke = h i [kh tM ^ (A. 14a) 

U l o jku = 0, (A. 14b) 

V? m = + 2K t[k K e]J , (A. 14c) 

V» kn = 2D {i K j]k . (A.14d) 

Inserting this result into Cf_, one obtains 

CI = -DiDjQ ij - 2/"' A ;; - Q ij K k K kj 

-Vh(R + K 2 - K t] K lj ) - (Q ij + 2Vhh ij )V£ ju . (A. 15) 

Clearly the variable V^- u plays the role of a Lagrange multiplier for the constraint + 
2\fhW = 0. The situation here is close to that discussed in Sec. 3e in which a constraint 
(denoted C\) is given by a momentum variable (denoted pi). In fact, the present theory can 
be placed in this form by a canonical transformation in which Q tj is replaced by Q tj + 2^/hh lj 
as the momentum conjugate to K^. The form of the canonical transformation can be 
deduced from the relationship 

/"'/',, + Q'-'K,, 

= (P ij + VhKh ij - 2VhK ij )h l3 + (Q l] + 2y/hti j )kij - 2(VhK)- . (A.16) 
Thus, define the new momenta 
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pij = pij + y/hKW - 2VhK ij , (A.17a) 
q« = gij + 2Vhh ij , (A. 17b) 

and the action becomes 

S E = iJ dt ^d 3 x{P ij h ij + Q l3 K l3 - 2{VhKy - NCf_ - V*Cf) 

+i J dt J d 2 x( v ^/Vh)n i {-2P ij Vj - 2Q ij K jk V k + Q ij DjN 

-NDjQ ij - 2y/hD i N + 2VhKV i ) , (A.18) 

where 

CI = I);!)/)" - 2I> : 'K U - Q ij K k K kj 

-Vh(R - K 2 + K i3 K lj ) - Q ij V° ju . (A.19a) 
Cf = -2D^' + Q jk DiK jk - 2Dj(Q jk K ki ) . (A. 19b) 

The variables Q 1,7 and (which play the role of pi and A 1 in the discussion of Sec. 3e) 
can be eliminated by setting Q l i equal to zero. The action then reduces to 

S E = i J dt J^d 3 x{P ij h ij - 2(VhK)' - V i 2i) ! i>' 

-N[-2P ij K i3 - - Vh(R -K 2 + KijRV)] } 
+i J dt J d 2 x^n { {^pt'Vj/Vh - 2D i N + 2KV 1 } , (A.20) 

where K i3 (which plays the role of q 1 ) appears as an extra independent variable in addition 
to the canonical pair h i3 , P tj , the lapse function N, and the shift vector V 1 . 
The equation of motion for K i3 which follows from the action (A.20) is 

= -2P ij + 2VhKh ij - 2VhK ij , (A.21) 

and the solution of this equation is 

y/hKij = -P^ + Phij/2 . (A.22) 

By substituting this result into Eq. (A.20), we find the action for Einstein gravity in the 
Hamiltonian form 

S E = iJ dt J^d 3 x{-P 3 h tj - V l \-2D 3 Pl\ - N[{2P lj P tj - P 2 )/{2Vh) - VhR] } 

+i J dt J d 2 x^n, [V 3 (Ph ij - 2P ij )/Vh - 2D*n} , (A.23) 

— ij — . . • 

Note that the "kinetic" term of Eq. (A.23) is —P h i3 , rather than the more usual P ZJ hij. 
The difference is just a boundary term, % / s d 3 x P at the initial and final times. If we had 
originally chosen the action (A. 9) to include a boundary term 2i f j: d 3 x VhK at the initial 
and final times, then we would have obtained P^h^ for the kinetic term in the Hamiltonian 
form of the action. 
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APPENDIX B: SURFACE GRAVITY 



As described in Sec. 4, one wedge of the spacetime g ab is foliated into stationary hypersur- 
faces t = const with time flow vector field t a , where t a is the Killing vector field that vanishes 
at the bifurcation surface TC. In this appendix I will drop the tildes with the understanding 
that all relationships hold for the Lorentzian metric g a b- 

The lapse function N and shift vector V a satisfy 

t a = Nu a + V a . (B.l) 

This expression is well defined on the interior of the wedge. Consider the limit in which TC is 
approached from within a t — const hypersurface, say, S . Since t a vanishes at TC, it follows 
by contraction of Eq. (B.l) successively with u a and h b a that 

limiV = 0, liml/ a = 0. (B.2) 

This result is used in Sec. 4 to show that the "terms ~ N and V in vanish at the boundary 
TC of E. 

The surface gravity of a Killing horizon equals |22[ 



k = Um(|t||a|) , (B.3) 

where \t\ = \/—t a t a is the magnitude of t a and \a\ — ^a c a c is the magnitude of the acceler- 
ation a c = (t a V a t c )/\t\ 2 of the orbits of t a . I will now show that, in effect, the shift vector 
V a vanishes sufficiently rapidly in the limit as the bifurcation surface is approached from 
within E so that the orbits of t a become orthogonal to the t = const surfaces. Then the 
surface gravity can be expressed in terms of the acceleration of the unit normal u a of the 
t = const surfaces. 

Using the Killing vector field property V (a ^) = 0, one can easily show that \t\a c = V c \t\. 
Then the surface gravity (B.3) becomes 



K = \imJ(V a \t\)(V a \t\) . (B.4) 



Observe that the limit of \t\ = \J—t a t a = ^/N 2 — V a V a is a constant (namely zero) as 7i 
is approached from within E . Thus we have limo" ab Vfe|t| = where <j a b = h a b — n a rib is 
the induced metric on TC and n a is the unit normal of Ti in E . Since also t a V a |t| = 0, the 
gradient of \t\ lies entirely in the v a direction, 

limV a |t| = \imv a v b V b \t\ , (B.5) 

where v a is the unit vector orthogonal to both a a b and t a : 

v a ~ lim(iVn a + n h V b u a ) . (B.6) 



Because the surface gravity of a bifurcate Killing horizon is nonzero [p2jl , it follows from 
Eq. (B.4) that lim V a |t| ^ 0. Also note that V a |t| is spacelike, since it is orthogonal to the 
timelike vector t a . Therefore lim /i^Vb|t| 7^ and, since limcr ab Vf,|t| = 0, we conclude that 
lim n b V b \t\ ^ 0. 
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Now consider the limit of V a / |t | as 7i is approached from within So- This is an inde- 
terminate form, 0/0. We can apply l'Hopital's rule and differentiate both numerator and 
denominator along the normal n a direction within S . The derivative of the denominator, 
n b V(,|t|, has a nonzero limit by the argument above. The derivative of the numerator is 



n b V b V a = n b V b {h a c t c ) 

= n b t c V b h a c +n b h a c V b t c . (B.7) 

Using h a c = 8% + u a u c , one can write the first term as 

n b t c V b h a c = n b t c V b {u a u c ) 

= NK ab n b - u a n h K bc V c , (B.8) 

where K ab = —h^V c u b is the extrinsic curvature of S . The two terms in Eq. (B.8) vanish 
in the limit because the lapse N and shift V c both vanish in the limit. Using the relation 
h ac = a ac + n a n c and the Killing vector field property V ( bt c) = 0, one can write the second 
term of Eq. (B.7) as 

n b h ac V b t c = n b a ac V b t c 
= -n b a ac V c t b 
= -n b a ac V c (Nu b + V h ) 

= Na ac K cb n b - n b a ac V c V b . (B.9) 

The first term in Eq. (B.9) vanishes in the limit due to the factor of N, and the second 
term in Eq. (B.9) is zero in the limit since the derivative acts along the bifurcation surface 
where V b vanishes. The result is that \imn b 'V b V a = 0, so by l'Hopital's rule we have 
limV7|t| = 0. Since \t\ = VN 2 - V b V b , we also find \im(V a /N)(l - V^/N 2 )- 1 ' 2 = 
which implies lim\^ a /iV = 0. 
The gradient of \t\ is given by 

V„| (H ^ + ^. (B.10) 

The results above show that the second term vanishes in the limit as TC is approached from 
within So, and that limiV/|t| = 1. Therefore Eq. (B.10) yields limV a |t| = limV a iV. We 
also find from Eq. (B.6) that v a = n a . Thus, Eq. (B.5) becomes 

lim V a iV = lim n a n b V b N . (B.ll) 

It follows that the surface gravity (B.4) can be written as k = lim|n a V a iV|. If we choose, 
as in the main body of the paper, the unit normal n a to point "radially inward" towards 
Tl, the surface gravity becomes k = — lim n a V a iV. From the relationship (B.ll) we find the 
key result 

hmV a iV = -hm/«n a . (B.12) 

Finally, note that the surface gravity also can be expressed as k = — lim Nn b A b , where (see 
Ref. f2T)|| ) A b = u a V a u b = hl(V c N)/N is the acceleration of the unit normal of the t = const 
hypersurfaces. 
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